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Abstract 

We investigate the mechanism of extraordinary optical transmission in subwave- 
length metal hole arrays. Experimental results for the arrays consisting of square or 
rectangle holes are well explained about the dependence of transmission strength on 
the polarization direction of the incident light. This polarization dependence occurs 
in each single-hole. For a hole array, there is in addition an interplay between the 
adjacent holes which is caused by the transverse magnetic field of surface plasmon 
polariton on the metal film surfaces. Based on the detailed study of a single-hole 
and two-hole structures, a simple method to calculate the total tranmissivity of hole 
arrays is proposed. 

PACS numbers: 78.68. +m, 78.20.-c ,42.25.Hz, 78.35. +c, 42.97.-c 
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1. Introduction 

The extraordinary optical transmission (EOT) [1,2] in a subwavelength metal hole array is 
an interesting topic [3- 12] because its mechanism is still in exploring and it shows abundant 
features. One of the features is that the transmissivity may depend on the polarization 
direction of the incident light. Disclosing clearly the reason behind the dependence is helpful 
to adjusting the EOT strength, as well as to applying the EOT in optical devices. Lots of 
experiments have been done for light in visible [3-7,13-17], infrared [8-10], and terahertz 
[11,12,18,19] regions to observe the dependence of EOT on polarization of the incident light. 
We here sort both the arrays and holes into three kinds, respectively, as summarized in 
Table 1. It is seen from Table 1 that among nine structures, five have been fabricated to 
observe the dependence of the EOT on the light polarization. In the second column, a square 
lattice consisting of rectangle holes shows polarization dependence while that consisting of 
square holes does not. In order to disclose the reason behind the discrepancy, a theoretical 
investigation is desirable. 



Table 1. The experimental results that whether EOT is dependent on the polar- 
ization of the incident light or not. The names in parentheses are used in Sec. 5. 





Square lattice 


Rectangle lattice 


Single hole 


Square hole 


Independent [11,16] (S-S array) 


Unreported (S-R array) 


Unreported 


Rectangle hole 


Dependentfll, 15-19] (R-S array) 


Unreported (R-R array) 


Dependent [14,15] 


Circle hole 


Indep endent [3- 5 , 1 3] 


Dependent [5-7,11,13] 


Unreported 



A lot of theoretical works about EOT in hole arrays have been reported[20], they are 
mainly focused on the mechanism or factors that cause in or influence EOT in hole arrays. 
A few of them investigate the polarization dependence of EOT in hole array or single hole 
[21-25] . Garcia et al. [21] carried out a rigorous solution of Maxwell's equations so as 
to obtain the transmission of circle holes perforated in a thin perfect-conductor screen for 
s and p polarization. Gordon et al [22] explained the polarization dependence in array 
of elliptic holes in terms of the interaction between SPP and the periodic lattice grating. 
Notwithstanding the approaches done, the systematical investigation about the polarization 
dependence is still desirable. 
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In our opinion, when considering the EOT in an array consisting of holes (or slits), there 
are basically the single-hole (slit) effect and the inter-hole (slit) effect [26]. The former 
reflects the transmissivity behavior of the light going through single subwavelength holes, 
and the latter means the possible modulation of transmissivity arising from the influence by 
neighboring holes. Supposing that the polarization dependence does exist, then one should 
know if the dependence is caused by the single-hole or inter-hole effect or both. 

In this paper we investigate the mechanism of the polarization dependence of the EOT in 
a hole array consisting of square or rectangle holes. Based on our simulation results by finite- 
difference time-domain (FDTD) method [27], we reveal the mechanism of the polarization 
dependence and explain the experimental results. Furthermore, we find that it is possible 
to get a simple way to calculate the transmissivity of the hole array, which may avoid the 
burdensome simulation work in hole arrays. 

This paper is arranged as follows. In section 2 the hole-array model is established. 
Before studying the EOT of the array, we study in detail the EOT of a single-hole and 
double-hole structures in sections 3 and 4, respectively, so as to clearly show the single-hole 
and double-hole effects. Then the EOT in a hole array is researched in section 5. In doing 
so, a simple method is proposed to calculate the transmissivity of the hole array. In section 
6 the simple method proposed in Sec. 5 is applied to some arrays. It appears that the 
application is satisfactory. Section 7 gives our summary. 

2. The array model 

Our model is sketched in Fig. 1. Rectangle holes drilled on a metal film form a two- 
dimensional lattice, consisting of rectangle cells, in xy plane. The lengths of two sides of each 
hole are a and b, and the lattice constants are A and B, respectively. A linearly polarized 
TM wave illuminates this structure along ^-direction. Before impinging the structure, the 
magnetic and electric components are Ho and Eq, respectively. In simulation we always use 

E = H = 0.31. (1) 

The angle between y axis and Eq direction is 9. Hereafter the light is termed as " 
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^-polarized" . From Fig. 1 the x and y components of electric and magnetic fields are 



The metal film with thickness h = 2 pom is made of silver. The incident wave length is 
Ao = 0.6/xm . The dielectric constant of silver vs. wavelength can be expressed as e^g = 
3.57-54.33A§+i(-0.083A +0.921Ag)[26,28]. Thus, as A = 0.6fim,e Ag = -15.989+^0.1491. 

When light goes through the array, the surface plasmon polariton (SPP) will be excited 
in every hole and the metal surfaces. Since it is a TM wave, after entering the holes, the 
electric field may have a ^-component, while the magnetic field does not. In each hole there 
is a strong power, denoted as P. In simulation, this power value P is measured by a monitor 
M placed at the exit of the hole labeled by " 0" . If the whole structure is removed, the power 
measured by this monitor at the same place is denoted as Pq [20, 22]. The transmissivity of 
this hole is defined as T = P/Pq. 

Our simulated result is the transmissivity T. Since T is simply linearly proportional to 
Poynting vector S, we will analyze the construction of S to explain the expression of T 
obtained by simulation. 

Since the array consists of holes, the light behavior in any one hole and the correlation 
between holes are essential in realizing the light behavior when light going through the 
whole array. Therefore, before studying the whole array, we explore the light behavior when 
it goes through only one hole and a two-hole structure. 

3. The single-hole structures 

Letting all the holes in the array except the one labeled by " 0" in Fig. 1 be closed, we set 
up a one-hole structure. We will study the cases where the hole is a square and a rectangle, 
respectively. 

Suppose that the amplitudes of the components of electromagnetic field in the hole are 
E^ ole , Ey° le , H^ ole and Hy° le , respectively. It is found from the simulated results that these 




(2b) 



(2a) 
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amplitudes as functions of angle 9 can be expressed by following way: 

(3a) 



E hole = E hole s{n Q 



E hole = E hole cog £ 
jjhole = -Hhole cos Q 

H^ ole = H^ e sin 9 



(3b) 



Equation (3) tells us that, when studying a single rectangle hole, one merely needs to measure 
the field components at an arbitrary polarization angle 9 so as to get the amplitudes-E^ 6 , 
E^°y e , iJo°f e and i?Q°' e in the hole by Eq.(3). Then the field components at any other angle 9 
can be easily calculated in terms of Eq.(3). 

It is noticed that the angular dependence of the left hand side of Eq. (3) is identical to 
that of Eq. (2). Based on this fact, we introduce a concept of SPP polarization excitation 
ratios (PERs) of each field component: the ratio of amplitude of field component in the hole 
to that in vacuum. They are denoted by rjEx, VEy, Vhx, and r]Hy respectively, and expressed 
as follows: 

„ rphole I rp rphole I rp 

VEx - & x /-box - &Q tX /&0 

„ rphole 1 rp rphole I rp 

VEy - Ey /b 0y - h Q y /t,Q 

„ rrhole I rr rrhole I rr 

rjHx - ri x /H 0x - H x / H 

„ rrhole I rr rrhole I rr 

VHy - ti y /H 0y - H Q y I H 

For example, t]e x is called the x-component SPP electric field PER. Since the SPP electro- 
magnetic field PERs are independent of angle, they are used to describe the basic property 
of the hole. 

As long as the field amplitudes in the hole are measured, the total curve of the transmis- 
sion power or transmissivity can be obtained. To explain the transmissivity, one needs to 
calculate Poynting vector S = E x H. 

In the hole, the averaged value in a time period of z-component of Poynting vector is 



(4a) 
(4b) 



ghole = ( E hole H hole _ fiholejjholej ^ ^ 
ghole = ( E hole H hole ^2 g + £ hole R hole ^2 gy^ ^ 

Let us first investigate the case of a square hole. The parameters are taken as a = b = 
0.2/im. The power measured in this hole is denoted as P hole and the transmissivity of the 
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hole is defined as To = P hole / Pq. The simulated transmissivity values are plotted in Fig. 2 
by crosses. Figure 2 shows that as 9 angle changes, To varies between 19.31 and 19.38. Since 
the variation scope is within calculation error, the transmissivity is regarded as unchanged, 
i.e., it is independent of the polarization. 

The calculated amplitude values of x and y components of electromagnetic fields of SPP 
in the hole are presented in Fig. 2 by square and circle symbols. The results are easily fitted 
by curves in terms of Eq. (3). Here, we have 

E hole = E hole = Lg8| (6a) 

#o,f = Kv e = 1-21. (6b) 

Equation 6 exhibits two features. One is that the dimensions of the field component ampli- 
tudes E^ e , Ey°y e , iTj^and H^'m the hole are larger than those outside the hole E Qx , 
E 0y , H 0x and H 0y respectively. This reflects the EOT character, i.e., the transmissivity T is 
greater than 1, as a subwavelength hole should have. The other is that the electromagnetic 
fields in a square hole behave as isotropic. This is an important feature of a square hole, 
which has been discovered by experiments, as shown in Table. 1. Later we will see that for 
a rectangular holes it is not so. 

For a square hole, substituting Eq. (6) into Eq. (5b), one obtains 

ohole jphole Tjhole lr> jphole jjhole lr> ffy\ 
b z = h 0,x H 0,y l l = h 0,y H 0,x I Z ' {') 

Here S^^is independent of angle 9, which is the reason why the transmissivity T is inde- 
pendent of angle 9 as shown in Fig.2. 
Equation (6) leads to 



VEx = VEy 
VHx = rj H y 



(8) 



In other words, along the two sides of the square hole, the SPP PERs are equal. Equation 
(8) is the physical reason that the transmission power is independent of the polarization 
direction in a square hole. 

Equation (8) manifests the ir/2 rotation symmetry of a square hole. It is probably 
that if the n/2 rotation symmetry is broken, Eq. (8) will not be valid. Consequently the 
transmission power should change with the polarization angle. 
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Next we investigate the case of a rectangular hole. Experiments showed the polarization 
dependence of transmission in rectangular hole, see Table 1. Let us see our simulated results. 

The simulation results of a rectangle hole are plotted in Fig. 3 where the two sides of 
the hole are a = 0.2/xm and b = 0.1/xm, respectively. The crosses in Fig. 3(a) show the 
transmissivity T . It changes with the polarization angle. The symbols in Fig. 3(b) denote 
the field amplitudes in the hole, which can be fitted with Eq. (3) but E^ e = 0.69 x KT 2 , 
Effi = 2.89, H^ e = 1.21 , and H%° y le = 0.53 x 10~ 2 . With these data we obtained: 
Vex = 0.022, r)Ey = 9.32, t]h x = 3.90, and rjHy = 0.017. That is to say, 



The SPP PERs in x and y directions for either electric or magnetic field differs from each 
other. This feature is different from that of a square hole. In each direction, when the 
electric field is strong, then the magnetic field is weak, or vice versa. 

The above discussions demonstrate that the SPP polarization excitation ratios are the 
key roles to exhibit the properties of the transmission with the polarization angle. 

As has been mentioned above, one can choose an arbitrary polarization angle 9 to get the 
amplitudes Eq ^, Eq ^, ifg°? e and H^ 6 . In this way, one saves a lot of time and workload 
to avoid measuring the whole angle region as shown by solid and open symbols in Figs. 2 
and 3. 

Since for a square hole, the SPP PERs in two axis directions are the same, while for a 
rectangle one it is not, when changing the ratio b/a, the SPP PERs should vary. We keep 
a = 0.2/xm, and change b from 0.1 to 0.3 fim. The simulated SPP PERs are displayed in 
Fig. 4 by symbols. The lines in Fig. 4 are just for guiding eyes. When the ratio b/a is 
small, say for b = 0.1/im, one of r] x and rj y is negligible compared to another. It seems 
that in this case the field components E^ ole and Hy° le are totally depressed. This implies 
that, the SPP wave is mainly polarized with the magnetic field along the longer side of the 
rectangle hole. In other word, with respect to the polarization properties, the rectangle 
hole is somehow equivalent to a slit when the ratio b/a of the hole is small. For convenience 
we refer to this kind of hole as slit-hole. The character of a slit-hole is that the electric 
or magnetic field component along one side direction is negligible compared to the other 
side. According to our simulation results, there exists a critical size for b beyond which the 




(9) 
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PERs r\ x and r] y aie comparable to each other. For instance, when a = 0.2/xm, it is found 
from Fig. 4 that the critical size for b is b c = 0.18/xm. When b > b c , the corresponding 
SPP PERs rise suddenly so the character of a slit-hole disappears. In other words, the 
SPP mode of electric field in x direction begins to excite. As the size of b approaches that 
of a, 7] x and ^become closer. At b = 0.2/im, the square hole, the two solid lines meet 
at the value r)E X = VEy = 6.06, and the two dashed lines meet at the value rjHx = VHy = 3.89. 

4. The double-hole structures 

Now we turn to investigate the inter- hole effect. For this purpose we set up a two- hole 
structure by closing all the holes in the array except the two labeled by "0" and " 1" , hereafter 
referred to as 0-1 structure. 

Before starting the investigation of the two-hole structure, let us briefly retrospect the 
EOT of a double-slit structure [20]. The SPP wave excited in one slit will interfere with 
that coming from the other slit. The interference varies with the inter-slit distance D. 
This interference is the so-called inter-slit effect. As a consequence, the total power passing 
through the structure oscillates with the inter-slit distance D. Each peak of the power curve 
corresponds to the in phase interference between the two slits. 

A double-hole structure resembles a double-slit structure in that there is an interference 
between SPPs excited in the two holes as the SPP waves travel along the metal film surfaces, 
and the interference varies with the inter-hole distance. Therefore, at appropriate inter-hole 
distances, the interference will generate largest transmission power. In simulation, we find 
that one of such distances is A = 0.4/xm when a x b = 0.2 x 0.2/im 2 and the polarization of 
the incident wave is 9 = 90°. 

The transmissivity measured in this structure is denoted as Toi. When fixing A = OAfim, 
the variation of T i as a function of the polarization angle 9 is displayed by the solid circles 
in Fig. 5(a). The data are well fitted with the equations 

T 01 (9) = T + AT 01 (9), (10) 

where 

AT O i(0) = 6.64 sin 2 9, (11) 

and T = 19.36. It should be noticed that T is just the transmissivity of single square 
hole structure, see Fig. 2. Equation (10) clearly demonstrates that the total transmissivity 
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comprise two parts coming from the single hole and inter-hole effects, respectively. The 
contribution from the inter- hole effect on the transmissivity is expressed by Eq. (11). When 
9 = 0, the inter- hole effect vanishes. 

Here we intend to disclose which factors the coefficient 6.64 of the interference term 
involves. To do so we analyze the behaviors of the SPPs traveling from one hole to the other 
along x direction. That is to say, we should evaluate the Poynting vecter S^ pp between 
the two holes. For this purpose we choose the mid-point between the centers of the two 
holes on the exit surface of the film as an observation point to measure S^ pp . Now we 
close the hole "1". Thus the system degenerates to a single-hole structure. In this case the 
SPPs at the observation point comes from the open hole. The observed 5*f pp reads S^ pp = 

fiSPF 'jjSpp _ fiSPF 'ffSPP Q ur simulated results reveals that the absolute values of Ef pp 

y z z y y 

and iff pp are negligible compared to Ef pp and Hy PF '• Thus we have S~! pp = —Ef pp Hy PP . 
Consequently, in the following, we merely take into account the contribution of Ef pp and 
Hy PP , omitting those of Ey PP and iff pp . The simulated amplitudes of E^ pp and iff PP at 
the observation point are displayed in Fig. 5(b) by solid and open circles, respectively. The 
solid curves in Fig. 5(b) are plotted with 



E$ pp = Eg pp sm9 
Hl PP = Hl PP sin9 



(12) 



where E$ pp = 1.12 and H§f p = 0.82. Comparing Eqs. (12) and (3), we find that the 

y ±A 0y I ±J -0y 



ratio Hy° le /Hy PP = Hq° /Hq PP is independent of angle 9, so we define this ratio as 701: 



I H y SPP /H^ = H SPP /H^, whenH^ ^0 
7oi = < (13) 
I 0, whenH^° le = 

This ratio is regarded as the conversion of the transverse magnetic filed on surface excited 
by the magnetic field in the hole. It is determined by the geometry and physical parameters 
of the single-hole structure. The subscript "01" in 701 means the SPP traveling from holes 
to 1. In the present case, 701 = 0.82/1.21 = 0.68. When H^ le = 0, the H^ pp can not be 
excited, so we define 701 = 0. In this case the inter-hole effect vanishes. 

Now we are ready to disclose which factors the coefficient of sin 2 9 term in Eq. (11) 
comprises. Let both holes in the structure open. The SPP wave from hole will arrive at 
hole 1, it will interfere with the SPP excited in hole 1. The same process will also occur at 



hole 0. In this case the total magnetic field at the exit of hole "0" includes two parts: One 
is if^ e contributed from "0" itself as an isolated hole, and the other is Hy PP contributed 
from "1". Thus the y-component of magnetic field becomes Hy° le + Hy PP . Consequently, 
the Poynting vector in z-direction is 

Qthole j^hole ^^jhole _|_ fjSPP^j jjjhole jjhole ghole _|_ j^hole jjSPP (14) 

Here S^ ole is just Eq. (5a), the Poynting vector without the contribution of Hy PP . 
Substituting Eqs. (3) , (12), and (13) into Eqs. (14) we get 

ohole ohole , A ohole /-i r \ 

where 

AS™: = 2^701 sin 2 9, (156) 

The Poynting vector of the 0-1 structure Sq° 1 1 comprises two parts: a term of a single-hole 
structure and a term reflecting inter-hole effect between the two holes. The latter in turn is 
related to its self. Therefore, although an angle factor sin 2 9 is separated, the coefficient 
of sin 2 # in Eq. (15b) should generally still contain functions of angle 9. Equation (15) is 
linearly proportional to Eq. (10). Therefore, we can reasonably rewrite Eq. (11) in the 
following form: 

ATo 1 (^) = To 1 (^)C , oi(^)sin 2 ^, (16) 

Coi(#)is the coupling coefficient that reveals the strength of inter-hole effect generated by 
total transmission power. Combining Eqs. (10) and (16) one achieves 

T 01 (9) =T + T 01 (9)C 01 (9) sin 2 9, (17) 

We emphasize that Eq. (17) is applicable to the double- hole structure consisting of identical 
rectangle holes. 

In the case of square holes, we have, from Eq. (11), T 01 (9)C i(9) = 6.64. Thus the two 
expression of the two factors T i(^) and Cq\{9) are easily solved. 

C O1 (0) = 6.64/ (T + 6.64 sin 2 9), (18) 

and 

T 01 (9) = T /[1-C 01 (9)sin 2 9}. (19) 
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Please note that since To = 19.36, much larger than the term 6.64 sin 2 9,Coi(9) in Eq. (18) 
is approximately a constant. Indeed, if we select 9 = tt/2, then Cqi{9) = 0.255, and the 
calculated T i(9) is plotted in Fig. 5(a) by dashed line. Apparently, this is a quite good 
approximation. 

It is worthy to point out that the inter-hole effect involves the contributions from SPP 
waves of both surfaces of the metal film. Equations (10), (17), and (19) have included the 
contribution from both surfaces. 

When the nearest neighbor (nn) hole is farther and at another azimuth angle, the in- 
terference between the two holes will vary. As an example to demonstrate this, we choose 
the two holes labeled by "0" and "5" in the array depicted in Fig 1, while other holes are 
closed, to discuss the next nearest neighbor (nnn) inter- hole effect. The array constants are 
A = B = OAfim and a x b = 0.2 x 0.2/im 2 . The simulated transmissivity as a function of 
the polarization angle are plotted in Fig. 6 by solid points. The data are well fitted by solid 
curve in terms of following expressions: 

T 05 = T + AT 05 , (20) 

AT 05 = T o5 Co 5 (0) sin 2 (£ + 45°), (21) 

and 

C 05 (9) = -2.55/[T - 2.55sin 2 (# + 45°)]. (22) 

Compared to the 0-1 structure, the 0-5 structure shows differentia in two ways. One is 
that the phase shift comes from the fact that the hole "5" is located at azimuth angle 45°. 
Correspondingly, there is the same phase shift in Hy PP compared to Eq. (12): Hy PP = 
Hq PP sin(# + 45°) . Therefore when 9 = 135°, there will be no propagation of Hy PP between 
the holes "0" and "5", i.e., the inter-hole effect vanishes at this angle. Indeed, from Eq. (21) 
the interference term is zero at this angle. The other is that the figure —2.55 in Eq. (22) is 
in the place of 6.64 in Eq. (18). 

The discussion about the interference in the two-hole structures above only concerns the 
azimuth. Another factor affecting the interference is the distance between the two holes. 
Let the distance between the two holes be r. Then the transmissivity oscillates with r. This 
oscillation is embodied in the value of C. Our simulation results show that for present square 
lattice of A = B = 0.4/xm, as r = A, T i(9)Coi(9) = 6.64, which just corresponds to the 
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interference in phase; and when r = \/2A, this distance makes the interference out of phase 
so that the transmission is suppressed, thus T^Cq^{9) = —2.55 is minus. From Eq. (21) we 
see that as 9 = 45° the inter-hole effect term is maximum, so that the transmissivity curve 
in Fig. 6 shows a valley. 

We have mentioned that in the 0-1 structure, calculated Tqi(9) using a constant 
C(9) = 0.255 approximates the exact results quite well, as shown in Fig. 5. Here we again 
set a constant C{9) = —0.152 to compute T 05 (6 I ) and the results are plotted in Fig. 6 by 
dashed line, which is almost identical to the solid line. 

5. The hole arrays 

To simulate the transmission of a hole array is quite difficult for a very large memory 
size is needed. However, the discussion about the two-hole structures in Sec. 4 prompts us 
that an array can be regarded as a combination of two-hole structures. Here we propose a 
simpler method to treat the hole array. 

As an example, we first consider a three-hole structure consisting of the open holes labeled 
by "1", "0" and "3" while other holes being closed in the array depicted in Fig. 1, referred to 
as 1-0-3 structure. In such a structure, if the transmissivity of the hole "0", denoted as T103, 
is measured, one has to consider inter- hole effects between hole "0" and its two neighbors . 
Thus a reasonable expression should be 

Tx 03 (9) = T + AT 01 + AT 03 = T + T 103 (9)C 0l sin 2 9 + T 103 (9)C 03 sin 2 (# + 180°) (23) 

Note that similar to the cases of 0-1 structure and 0-3 structures, the coefficients of the two 
interference terms should include a factor of the total transmissivity T 10 3(#). Since the holes 
"1" and "3" are symmetric with respect to the hole "0", we have Cqi = Cq 3 . The difference 
of Tio3(#) between results of the calculation with Eq. (23) and simulation by FDTD method 
is about 1%. 

From the example of the 1-0-3 structure it is concluded that for each additional nn hole, 
one merely simply add a term to embody the inter-hole effect, although the coefficient should 
be proportional to the total transmissivity of the hole "0". This conclusion can be extended 
into the whole array. 

Now let the all holes in the array open. We calculate the transimissivity of hole "0". There 
are four nn and four nnn neighbors around this hole, labeled by " 1" to "8" , respectively. The 
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influence of the holes farther than the nnn ones is merged into the inter-hole effect between 
hole " 0" and the eight neighbors, so that it needs not to consider. Thus, the transimissivity 
of hole "0" reads 

8 

Tarray = Tq + AToj. (24) 

i=l 

The second term in Eq. (24) includes the contribution from its all 8 neighboring holes. 
When the hole array composes a square lattices with A = B. 

First we study the case where all holes in the lattice are square, referred to as S-S array. 
The transmission in hole "0" is denoted as Tss- 

The 0-1 and 0-5 structures have been studied in detail in Sec. 4. According to the con- 
clusions of the two structures, we easily put down the terms of inter-hole effects contributed 
from all the eight neighboring holes as follows: 

AT 01 = TssCoi sin 2 6, (25a) 

AT 02 = T S -sC 02 sin 2 (£ - 90°) = T S . S C 02 cos 2 9, (256) 

AT 03 = T S -sC 03 sm 2 (6 - 180°) = T S _ S C 03 sin 2 9, (25c) 

AT 04 = T S -sCo4 sin 2 (# - 270°) = T S - S C m cos 2 9, (25d) 

AT 05 = T S _ S C 05 sin 2 (# + 45°), (26a) 

AT 06 = T 5 „ 5 C 06 sin 2 (# - 45°), (266) 

AT 07 = T S -sC 07 sin 2 (# - 135°), (26c) 

AT 08 = T 5 „ 5 C 08 sin 2 (# - 225°). (26d) 

Here the azimuth of each neighboring hole is taken into account. Since the four holes "1" 
to "4" have the same distance away from "0", and the other four on the vertexes do so too, 
one naturally gets: 

Coi = C02 = C03 = Co4 ; (27a) 
Cos — C*06 = C07 — Cos- (276) 
Inserting Eqs. (25)-(31) into (24), we obtain 

T s _ s = T + 2T S _ S (C 01 + C 05 ) (28) 
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Obviously, the transmission is independent of 9 angle. This explains the experimental result 
of S-S array listed in Table 1. 

If the parameters of the holes and lattice are the same as those in Sec. 4, we have 
C 01 + C 05 = 0.103, T = 19.35, thus T s _ s = 24.37. 

Next we study the case where the holes are square and the lattice is rectangular, referred 
to as S-R array. The transmission in hole "0" is denoted as Ts-r. Since in this case B ^ A, 
we define an angle a: 

a = arctan(5/A). (29) 

The angular dependences of AT 0i , i = 1, 2, 3, 4 are the same as those in Eq. (25). One merely 
need to replace T 5 _ 5 in Eq. (25) by T S -r, to get the expression of AT 0i , i — 1, 2, 3, 4. However, 
since B ^ A, Eq. (27a) is not valid any more. We have following relationship: 

Coi — C*03 7^ C*02 = C*04- (30) 

As for the neighbors "5" to "8", the angular dependence of interference terms are written 
as 

AT 05 = T S - R C 05 sin 2 (9 + a), (31a) 

AT 06 = T S - R C 06 sin 2 (9 - a), (316) 

AT 07 = T s _ R C 07 sin 2 [9 - (180° - a)] = T S _ R C 07 sin 2 (9 + a), (31c) 

AT 08 = T 8 - R C W sin 2 [9 - (180° + a)] = T S ^ R C 08 sin 2 (9 - a). (31d) 

The distance between the hole "0" and any one of these four holes is the same as others. 
Hence Eq. (27b) is still valid. Inserting Eqs. (30) and (31) into Ts~r = T + Yli=i AT j, we 
have 

T S _R = T + 2T s ^r{[Coi sin 2 9 + C 02 cos 2 9] + Co5[sin 2 (£ - a) + sin 2 (£ + a)}}. (32) 

It is seen that the transmissivity of the S-T array depends on the polarization angle 9. 

Thirdly we discuss the case where the holes are rectangular, a x b = 0.2 x 0.1/im 2 and 
the lattice is square, A = B = 0.4/xm, the so-called R-S array. Since the transmission 
in each rectangular hole is dependent on 9, as manifested in Fig. 3(a), the single-hole 
effect is enough to cause the dependence of transmissivity Ts^^on the polarization angle. 
Besides, the inter- hole effect also influences T R _$- Apparently, in this case we again have 
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relationship Eq. (30). In the case of a x b = 0.2 x 0.1/im 2 , the hole is regarded as a 
slit-hole, as having been mentioned in the last paragraph of Sec. 3 in explaining Fig. 4. 
Considering Eqs. (15) and (13), for a slit-hole, we have ASq° 1 ^ = 2Sq° 1 ^Jqi sin 2 9 = 0, so 
AT i(9) = Tqi{6)Cqi{6) sin 2 6 = 0. Consequently Cqi = 0. Following the method as above, 
the expression of the transmission is obtained as follows: 

T R _ S = T Q {6) + 2T R . S [C 02 (9) cos 2 9 + 2C 05 (9) sin 2 45° cos 2 9], 

By our simulation, it is approximately that C 02 (9) = 0.07/ cos 2 9, C 05 (9) = 0.12/ cos 2 9. 
Thus we get Tr_s — 39.2 cos 2 9. The feature of T R _s curve is the same as the experimental 
result [17]. This result confirms the validity of our calculation method. 

Finally, for the case of a rectangular lattice comprising rectangle holes, R-R array, we 
can use the same method to discuss the transmission T R _ R . But we do not put down the 
formula. A qualitative conclusion is obvious. Since both the single holes and lattice are 
rectangular, it is sure that T R _ R depends on the polarization angle. 

6. Applications 

Up to now we have discussed the six cases in Table 1. The two kinds of single holes 
are studied in detail in Sec. 3 and the four kinds of arrays are investigated in Sec. 5. 
The mechanism of the polarization dependence of the transmission of each case is explicitly 
disclosed. A simple method is proposed to evaluate the transmissivity of the arrays. The 
physical meaning of this method is that it exhibits the transmission is mainly from two 
parts: the single-hole and inter-hole effects. The obvious advantage of this method is that 
it reduces the workload greatly compared to the simulation of the whole array. 

Among four kinds of arrays in Table 1, two, S-S and R-S arrays, have been investigated 
experimentally, while the other two, S-R and R-R arrays, have not. Here we employ our 
method to calculate the transmissivity Ts-r and Tr_# of S-R and R-R arrays. The numerical 
results are provided for someone to test. 

For an S-R array, we chose A = 0.4/xm, B = 0.3/im, and a x b = 0.2 x 0.2/im 2 . In 
this structure, Cqi = 0.255 is unchanged, and C02 and C05 have to be estimated. The 
transmissivities of 0-2 and 0-5 structures are simulated and fitted by T02 = T0 + O.I2T02 cos 2 9 
and T 05 = T - 0.08T 05 sin 2 [# - arctan(0.3/0.4)], so we get C 02 = 0.12 and C 05 = -0.08. 
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Then the transmissivity is expressed by 

T 5 _ R = 19.36/{1 - 2(0.255 sin 2 9 + 0.12 sin 2 9) + 2 x 0.08[sin 2 (# - 36.9°) + sm 2 (9 + 36.9°)]} 

The calculated curve is displayed in Fig. 7(a). 

For an R-R array, we take A = 0.4/im, B = 0.3/xm and a = 0.2fim, b = 0.1/xm. For such a 
slit-hole, Cq\ = 0. In 0-2 and 0-5 structures, the simulation results are T 02 = T (9) + 6 cos 2 9 
and T 5 = T Q (9) + 6.1 sin 2 (36.9°) cos 2 9 respectively. Here T (9) is the single hole transmis- 
sivity expressed in Fig. 3(a) with the formula T (9) = 24.3 cos 2 9. When writing formulas 
in the coupling form, T 02 = T Q (9) + Co 2 T 02 cos 2 9 and T 05 = T (9) + C^To^sin 2 (36.9°) cos 2 6*, 
we obtain C 02 {9) = 0.2/ cos 2 9 and C 05 {9) = 0.023/ cos 2 9. It is found that C 02 
andCos vary with 9 and cannot be regarded as constants now. This arises from 
that T (9) varies with 9 in this structure. Equation (27b) still holds in this lat- 
tice. Cq§{9) = 0.023/ cos 2 9 = Cw{9) = Cos{9). Thus the coupling equation is 
T r _r = T Q {9) + 2T r _ r [Cq 2 {9) cos 2 9 + 2C 05 (8) sin 2 (36.9°) cos 2 9}. The calculated T R _ R is 
plotted in Fig. 7(b). Comparing Fig. 7(a) and (b), we see that the variation scope oiT R _ R 
is larger than that of Ts-r, since the R-R array has a stronger anisotropy than S-R array. 
Especially, Tr^r can be zero at 9 = 90°. 

7. Summary 

We have investigated the polarization dependences of the transmission in square and 
rectangular lattices consisting of different subwavelength holes. The filed components and 
transmissivities of single-hole and double-hole structures are computed by use of FDTD 
simulation method. The behaviors of the transmission are explored and the corresponding 
mechanisms are disclosed. Our basic point of view is that the total transmissivity of a hole 
array is determined by the two basic factors: the single-hole effect and inter-hole effects. 
Based on the results of these structure, a compact method is suggested and applied to 
investigate the hole arrays. Our conclusions are summarized as follows. (1) The SPP PERs 
are key roles in single hole. In a square hole the SPP PERs along the two sides of the hole are 
equal, which leads to two consequences. One is that the SPP wave in the square hole is along 
the polarization direction of the incident light, and the other is that the amplitude of the SPP 
wave is in proportion to that of the incident light at any polarization angle. By contrast, in 
a rectangle hole, the SPP PERs is not isotropic, which results in that the amplitude of SPP 
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in a rectangle hole cannot reserve a fixed proportion to that of incident light with different 
polarization angle. Therefore the transmissivity depends on the polarization angle. (2) The 
transverse magnetic filed of the SPP wave on the metal film surface plays a key role in the 
inter-hole effect. (3) The total transmissivity of the hole array can be expressed with the 
single-hole transmissivity plus the terms reflecting inter-hole effects between the hole and 
its nn and nnn neighbors. (4) The conclusion (3) provides a simple method to calculate the 
transmissivity of the hole arrays. By this method we calculated the polarization dependence 
of the transmissivity for S-R and S-S arrays that are not reported in literatures. 
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Figure captions 

Fig.l The sketch of a metal hole array consisting of subwavelength rectangle holes. 

Fig. 2 The simulated transmissivity T and polarization dependence of amplitudes of the 
electric and magnetic fields in a single-square-hole structure. The fitting curves are 
calculated by Eq. (3). 

Fig. 3 The polarization dependence of quantities for a rectangle hole, (a) The transmissivity. 
The fitting curve is Tq(8) = 24.3 cos 2 9. (b) The amplitudes of the electric and magnetic 
fields. The fitting curves are calculated by Eq. (3). 

Fig. 4 The SPP polarization excitation rations of a rectangle hole with a = 0.2/im and b 
varies from 0.1 to 0.3/im. The lines are just to guide eyes. 
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Fig.5 (a) The polarization dependence of the transmissivity T m for a two-hole structure. 
The fitting solid curve is from Eqs. (10) and (11). The dashed curve is obtain by the 
approximation of Coi = 0.255. (b) The amplitudes of the electric and magnetic fields 
along the metal surface yielded from a single hole. The fitting curves are from Eq. 
(12). 

Fig.6 The simulated transmissivity in a 0-5 structure. The dashed curve is plotted by the 
fitted formula in the approximation of C05 = —0.152. 

Fig. 7 The calculated transmissivity in (a) the S-R array and (b) The R-R array. 
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